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We introduce a new slave-boson mean-field theory 
which allows the investigation of general multi-band 
Hubbard models. Unlike earlier attempts of such a gener- 
alisation, in our approach the quantum-mechanical prob- 
lem is exactly reformulated in an extended Hilbert space 



of Fermions and Bosons before a mean-field approxima- 
tion is applied. Systems with superconducting order pa- 
rameters are naturally included in our formalism. Our 
ground-state energy functional agrees with the corre- 
sponding quantity derived within the Gutzwiller theory. 
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1 Introduction For the investigation of correlated elec- 
tron systems in two or three dimensions only few theo- 
retical approaches exist which lead to sensible results for 
small, medium, as well as strong Coulomb interaction pa- 
rameters. A widely used approach which fulfils this crite- 
rion is the slave-boson mean-field theory. It has been intro- 
duced by Kotliar and Ruckenstein for an investigation of 
the single band Hubbard model | l The approach is based 
on two main ideas: 

i) The original quantum-mechanical problem, described 
by a fermionic Hamiltonian H, is reformulated by in- 
troducing bosonic degrees of freedom at each lattice 
site i of the system. As long as certain constraints are 
exactly fulfilled, the resulting new Hamiltonian H in 
the enlarged bosonic-fermionic Hilbert space is mathe- 
matically equivalent to H. 

ii) The Hamiltonian H is investigated by means of a mean 
field theory or, in the language of functional integrals, 
by a saddle-point approximation. 

Unfortunately, there is an infinite number of Hamiltonians 
H , equivalent to H, which may all lead to different results 
on mean-field level. Therefore, the approach requires a so- 
phisticated guess in order to find the 'right' Hamiltonian 
which yields sensible mean-field results. The choice, made 
by Kotliar and Ruckenstein, leads to the same ground-state 
energy functional as it is found by an exact evaluation of 
Gutzwiller wave-functions in the limit of infinite spatial di- 
mensions Q3). 



For the investigation of real materials one usually has 
to take into account the multi-orbital electronic structure of 
such systems. This requires the study of multi-band Hub- 
bard models. While the generalisation of the Gutzwiller 
theory for multi-band models is rather straightforward even 
for systems with superconducting ground states |4][5][6]], 
the very same generalisation of the slave-boson theory turned 
out to be difficult. Only recently such a generalised slave- 
boson scheme has been introduced for the treatment of 
general multi-band models (7). The ground-state energy 
functional, derived in that work, has been shown to agree 
with the Gutzwiller functional [8 |. In a more recent work 
|[9l , the scheme introduced in |!7| was generalised in or- 
der to study superconducting systems. For such systems, 
however, the ground-state energy functional derived in [|9| 
seems not to agree with the corresponding Gutzwiller func- 
tional 0. 

In the derivations of Refs. QUI there is a fundamental 
shortcoming. While in the single-band case the derivation 
of the slave-boson theory is exact up to the point where 
the mean-field approximation is applied, the multi-band 
derivation of Ref. (7J requires additional approximations 
already on the level of the operator equations. It is the 
main purpose of this work to introduce an alternative slave- 
boson scheme for multi-band models which avoids these 
additional approximations and is exact apart from the fi- 
nal mean-field treatment. In addition, our new approach 
automatically covers systems with superconducting order 
parameters. 
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The presentation is organised as follows: In section [2 
we introduce the general multi-band models that are inves- 
tigated in this work. We briefly summarise the results of the 
Gutzwiller theory for multi-band Hubbard models in sec- 
tion [3] A reminder of the Kotliar-Ruckenstein theory for 
the one-band Hubbard model is given in section 14.11 Pre- 
vious attempts to generalise this approach are discussed in 
section 14.21 In section [5] we derive our new slave-boson 
theory for the investigation of general multi-band models. 
A summary closes our presentation in section [6] 

2 Model Hamiltonians We study the general class of 
multi-band Hubbard models 

" EE**/<«A>' ' E "<•'•■•• 

i^j a, a' i 

where we introduced the combined, spin-orbital indices a, 
hopping parameters t^ ,<T , and local Hamiltonians 



,ct c. 



(2) 



E * 



i,<7l 2,(T2 i,&3 i,(T4 



(71 , (T2 ,0"3 ,CT4 



for each lattice site i. The Hamiltonian (f2]) is determined 
by the orbital-dependent on-site energies £i; ffl <7 2 and by 
the two-particle Coulomb interaction jj°' 1 <' T2 >' 73 '' 74 _ it con- 
tains 2A^ local spin-orbit states cr, which we assume to be 
ordered in some arbitrary way, a = 1 , . . . , 2N. Here, N 
is the number of orbitals per lattice site. In order to set up 
a proper basis of the local Hilbert space, we introduce the 
following notations for the 2 2N possible configurations: 

i) An atomic configuration / is characterised by the elec- 
tron occupation of the orbitals, 

Ie {0;(1),...,(2JV);(1,2),...,(2,3), (3) 
...(2N - 1,2N);... ;(!,..., 2N)}, 

where the elements in each set / = (<n, 02, • • •) are or- 
dered, i.e., it is <7i < 02 < • • ■■ The symbol in (01 means 
that the site is empty. In general, we interpret the indices 
/ as sets in the usual mathematical sense. For example, 
in the atomic configuration I\I' only those orbitals in / 
that are not in I' are occupied. The complement of / is 
I = (1,2,..., 2N)\I, i.e., in the atomic configuration I 
all orbitals but those in / are occupied. 

ii) The absolute value |/| of a configuration is the number 
of elements in it, i.e., 



|0| = O;|(<7i)| 

(ci,er 2 )| = 2; 



1: 



(4) 



\(1,...,2N)\ 



2N 



iii) A state with a specific configuration I is given as 

\i) = c\ |o> = n4lo> = 4, •••4, jo) , (5) 



where the operators & a are in ascending order, i.e., it is 
<j\ < a-2 ■ ■ ■ < cr\i\- Products of annihilation operators, 
such as 

C I = Y[c a = c ai ...c asij , (6) 

will always be placed in descending order, i.e., with <ti > 
(72 ■ . ■ > <7|/|. Note that we have introduced the operators 

Cj and Cj just as convenient abbreviations. They must not 
be misinterpreted as fermionic creation or annihilation op- 
erators. 

iv) The operator rhj j> = \I) describes the transfer be- 
tween configurations I' and /. It can be written as 



miv 



c\c v n (1 



(7) 



cr"e.J 



where J = I U P. A special case, which derives from ©, 
is the occupation operator 



mi = \I) (I\ = l[nrl[(l 

a£l cr'el 



(8) 



The states \I) form a basis of the atomic Hilbert space. 
Therefore, any other basis \T) of the atomic Hilbert space 
can be written as 



\r)=J2Ti,r\I) 



(9) 



with coefficients Tj p. With such a general basis, the atomic 
Hamiltonian has the form 



r,r> 



m t . r .r> 



(10a) 



(10b) 



In case that we deal with only one orbital per lattice 
site, the Hamiltonian (Q]i reads 



Hib = E E ^4,<A<r + E U l rh l , 12 . (11) 

i 7 j a—1 i 



Here, the indices a = 1, 2 represent the two possible spin 
directions and rhi : i2 = rc*i,iHj,2 with hi, a = cj c i a is the 
'double-occupancy operator' on lattice site i. 



3 The Gutzwiller Energy Functional for Multi-Band 
Systems For later comparison with the slave-boson mean- 
field theories, we briefly summarise the results of the Gutz- 
willer theory for multi-band Hubbard models. For all tech- 
nical details, we refer the reader to Refs. |H1|5|. 

Multi-band Gutzwiller wave-function have the form 



i*G)=Aii*b>=n^i* r o> 



(12) 
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where \Wq) is a normalized single-particle product state and 
the local Gutzwiller correlator is defined as 



(13) 



Note that, instead of $13[ , one can also work with a non- 



diagonal variational parameter matrix A^ r ,, 

P i = J2 x r,r'\r)ii(r'\. 



(14) 



r,r' 



Since we work with an arbitrary atomic basis \T)i, how- 
ever, both correlation operators ( fT3l and ( TBI define the 
same variational space. In the following, we summarise the 
main results for a correlation operator of the form ( TP3V 
In general, the local uncorrected density matrix 



(c[<Aa'}#0 



(15) 



is non-diagonal with respect to cr, a' . In superconducting 
systems, it additionally exhibits anomalous elements such 
as (cj cj a }y? . By means of a unitary transformation (i.e., 
a Bogoliubov transformation for superconductors) one al- 
ways finds a local basis with a diagonal density matrix and 
vanishing anomalous elements. In the following, we only 
work with such a local basis, since it simplifies the results 
for the variational energy. Nonetheless, we use the nota- 
tions which we introduced in the previous section. For su- 
perconducting systems, this means that the single-particle 
Hamiltonian in (Q~|i has the more complicated form 



H 



EE 



[t 1;CT,CT 'ct c +t 2 '' a ' a 'c c f (16) 



,3;cr,cr' -t ~t i j4;<t,(t' ~ 

t C C -r- t C C 



Note that for our 'orbital' basis \a) the expectation value 
of the operator (|7]i with respect to |*£b) is given as 

(17a) 
(17b) 



where 



rr£l ail 



(17c) 



The variational parameters X r need to obey certain 
constraints, which naturally arise in the evaluation in in- 
finite dimensions HO. These are 

1 = ^K,rMr(™>i;r)v , (18a) 
r 

(c{<T%o->)*o = ^Kr^-siclaC^rhi-r)^ , (18b) 
r 

With the expectation value 

m i; r,r> = {w>i;r,r')& G = (mi-r,r')^ K,r^r' (19) 



where 



;r,r')f - y^Ti;j,r^.j,r' m i;/ 



(20) 



one can readily calculate the local energy (Hi.\ oc )$, G . F° r 
the calculation of the single-particle energy, we have to de- 
termine the expectation values of normal and anomalous 
hopping operators, 

(ee 2 ), G -E(^) w (^)*^5 ; ^)^' 

(7' lt a' 2 

(21) 

where, to simplify the notation, we dropped the lattice site 
index of the 'renormalisation matrix' q° . This renormali- 
sation matrix is given as 

<£' = ]T \* r \ r ,(r\cl \r')J2Ti, r Tp, AHfrha , 



r,r> 



where we introduced the operator 



H 



(l-/ CTV )(/ / |cjru<7 , >m /w 



(22) 
(23) 



+ (A^'|c CT ,|/) [U'.rm^j, + (1 - f a >j>)ml x<J , }1 



Here, we use the abbreviation f a j = (I\c^c a \I) and the 
operator 



<r'e J\<r 



(24) 



which is defined for a 6 J = I L) I'. Note that the ex- 
pectation value of d23l l in Eq. (1221 can be readily evaluated 
with equations ( TP7I ). 

4 The Kotliar-Ruckenstein Theory In the first part 
of this section, we introduce the auxiliary particle method, 
which was proposed by Kotliar and Ruckenstein for an in- 
vestigation of the single-band Hubbard model. In the sec- 
ond part, previous attempts to generalise this approach for 
multi-band Hubbard models are discussed. 

4.1 The One-Band Model We start from the Hamil- 
tonian ( fTTT i with its four-dimensional local Hilbert space 
Hi for each lattice site i, represented by the four states 
|J) = |0), \a) ,|12) (with a = 1,2 for the two spin di- 
rections). The Hilbert space of the whole lattice system is 
given by the tensor product 



H 



(25) 



Kotliar and Ruckenstein introduced auxiliary bosonic 
operators <jy. T , ^.j, which lead to an enlarged local Hilbert 
space Hf B defined by the basis states 



I,I%FB = \I)i»\I') 



r.Ii 



(26a) 



Copyright line will be provided by the publisher 



4 



J. Bunemann: Slave-boson theory for multi-band models 



Here, |i)jis the fermionic configuration state, defined in 
equation qJ, and |7) ,. B is the bosonic state 



Was = 4>lj |o>< 



(26b) 



with the bosonic vacuum state |0).. B . The original quan- 
tum mechanical problem can be recovered in the following 

way: 

i) One has to find a subspace Hi of Hf B which is isomor- 
phic to the physical Hilbert space Hi for each lattice site i. 
Kotliar and Ruckenstein defined this subspace by means of 
the constraints 



F, 



i.O 



1 



o, 



H:12 



with the bosonic occupation operators 



(27a) 
(27b) 

(28) 



The constraints (l27l i define the subspace Hi via the condi- 
tions 

Kb \t) = (29) 

for each G "Hi and <j G (0, 1, 2). Alternatively, we can 
define Hi directly by specifying its basis 



FB 



(30) 



The corresponding Hilbert space for the lattice system is 
given by 

H = <E)Hi . (31) 

i 

Note that, by construction, there is now a one-to-one cor- 
respondence of all physical states £ H and their coun- 
terparts \}P) G H. 

ii) With the auxiliary Hilbert spaces Hi and H properly de- 
fined, one can find operators Qi in H% that are similar to the 
physical operators Oi in Hi- Here, 'similarity' means that 



'Mm = iimin 



(32) 



for all configurations \I), \I'). With similar local operators 
Qi, one can set up an 'effective' Hamiltonian Hib which 
is similar to the physical Hamiltonian Hib, i.e., it obeys 



(V\H 1B \V') = {t\MiB\t') 



(33) 



for all physical states , G H and their counterparts 
W) i W) G %■ In this way, we nave introduced an exact 
mapping of the original physical problem, described by the 
Hamiltonian H\b in its Hilbert space H and the effective 
Hamiltonian Hib in H . 

To set up iJiB, we start with an identification of oper- 
ators that are similar to the fermionic operators c^l in Hi. 
Their counterparts cfl in Hi can be chosen as 



(34a) 



where the bosonic operators 

A —It 1 , 11 . 



i,o — Yi:12^i;a ' Vi;uVi-$ 1 

t — 1\ X L It 1 



z,<7 ri:<7 



J i;12 



(34b) 
(34c) 



have been introduced. As required, the operators cj^ obey 

equation ([32}. To set up the Hamiltonian Hib in Hi, we 
further need to find an operator m^i2 that is similar to 
rhi-,12 = hi^fii,2- The most obvious choice is 



TOi;12 = n®A2 ■ 



(35) 



Note, however, that there is a large amount of arbitrariness. 
For example, the operators 

m-i-12 = nf.jThi.i2 Or TO,; ; i2 = TOi;12 (36) 

are also similar to m.; ; i2 since both obey equation (l32b . The 
same ambiguity arises for the operators ( I34al ). For exam- 
ple, they may equally well be chosen as 



with 



-,.n 'L.rr' , 

=(l-A, CT )- 1/ M,.(4.)- 1/2 



and 



(37a) 

(37b) 
(37c) 

(37d) 



In fact, this choice is better than ( I34ab and was used by 
Kotliar and Ruckenstein since it yields the correct ground- 
state energy in the uncorrelated limit U = if the resulting 
effective Hamiltonian 



En 



^AiAA/ ^ + ^X>?i2 (38) 



is investigated on a mean-field level, see below. 

Kotliar and Ruckenstein | fl used a functional integral 
approach to calculate the free energy of the Hamiltonian 
(13 8t . For ground-state properties, i.e., at zero temperature, 
their saddle-point approach is equivalent to a replacement 
of the bosonic operators 7 by the amplitudes ipi-j. They 
govern the bosonic occupations 



H-I 



(39) 



and have to be determined by a minimisation of the ground- 
state energy functional 



(40) 



Here, the factors q* a and q i CT are defined in d37l i with the 

operators tpVj replaced by ip^.j. Note that all quantities in 
this section are real and the asterisks, e.g., in equation i40i . 
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are only used in anticipation of the corresponding multi- 
band results in section 

Instead of dealing with the exact constraints ( 1271 ). they 
are also satisfied only on a 'mean-field level' by Kotliar 
and Ruckenstein, i.e., for the expectation values 

1 = X>* J; (41a) 

< CT = ^12 + ■ ( 41b > 

The constraints (HTI) and the energy functional ( f4Qb are 
the same as those derived for the Gutzwiller wave function 
in the limit of infinite spatial dimensions II 1 Oil 111 or eval- 
uated by means of the Gutzwiller approximation 11211131 

m. 

4.2 Multi-Band Hubbard Models A generalisation 
of the slave-boson theory is straightforward 115111611 for 
multi-band Hubbard models with a local Coulomb inter- 
action of the form 

Hi = 22 U^ a ih a h a i = 2J Ujttij , (42) 

a, a' I 

where 

Ui= J2 U °>°' (43) 

a,u'£l 

It yields the same energy functional as derived within the 
Gutzwiller theory ifTTl . For the treatment of general multi- 
band Hubbard models a generalised slave-boson theory has 
been derived by Dai et al. ifTHl and, more successfully, by 
Lechermannet al. 0. 

As demonstrated in the previous sections, the slave- 
boson approach contains a number of adjustable objects. 
These rather flexible elements of the theory are the defini- 
tion of the extended Hilbert space Hf B , the definition of its 
physical subspace Hi, the form of the constraint equations 
Fie = and, finally, the particular definition of similar 
operators Qi. Despite this huge flexibility, in both works 
ifTHl and (7), the authors fail to derive an exact mapping of 
Hilbert spaces and Hamiltonians which, on a mean-field 
level, leads to satisfactory results. We will briefly sum- 
marise the previous attempts to formulate a generalised 
slave-boson scheme for multi-band Hubbard models in this 
section. Our own derivation for such an approach is dis- 
cussed in section [5] 

Dai et al. use constraint equations which do not de- 
fine the correct physical Hilbert space. This problem has 
been pointed out and solved by Lechermann et al. In their 
work, however, they fail to derive proper fermionic oper- 
ators 6^ . Instead, symmetry arguments are used in order 

to guess the form of certain operators cp , which lead to 
a reasonable energy functional on mean-field level. These 
operators, however, are not similar to the physical opera- 
tors c| , see below. Therefore, the whole derivation seems 
even less controlled than for the single-band model. In this 



section, we briefly summarise the main ideas of the slave- 
boson mean-field theory introduced by Lechermann et al. 

As in case of the single-band model one has to set up 
a local Hamiltonian Hi which is isomorphic to the physi- 
cal fermionic Hamiltonian Hi with its basis | r).. Lecher- 
mann et al. discuss various possibilities to define such local 
Hilbert spaces consisting of bosons and fermions. At first 
sight, it seems that Hi is most naturally defined as a gener- 
alisation of (130b through a basis 

= \n i ® \r)i-,B (44) 

with bosonic operators (jy^ r and the corresponding states 

\r) i . B = <£t r |0) i . B . However, the Hilbert-space Hi de- 
fined by this basis is discarded by Lechermann et al. on 
the grounds that there is no way to find a reasonable set 
of constraint equations as an alternative definition of Hi- 
In section |3J we show that the basis d44l > can, in fact, be 
used for a slave-boson theory, which, however, has to be 
different from the original scheme introduced by Kotliar 
and Ruckenstein. 

Instead of d44l i. Lechermann et al. introduce the basis 

\ r )i = -m £ 4Wio)i; B ®i^ (45) 

Vr I /(|/|=|r|) 

as a definition of their Hilbert space Hi- Note that in their 
derivation they draw a distinction between physical parti- 
cles and quasi-particles, described by operators dfl and 

fiJ., respectively. Our derivation in this section indicates 
that this distinction is unnecessary. 

The Ansatz (|45T > employs bosonic operators <j>i : rj for 
each pair of multiplet states \T) and configurations states 
| J) with the same particle number \T\ = \I\. The num- 
ber of these operators is much larger than the dimension 
of the local Hilbert space, which is different from the orig- 
inal single-band scheme introduced by Kotliar and Ruck- 
enstein. Despite this large number of operators, the space 
Hi is isomorphic to the physical Hilbert space Hi with its 
basis \r)-. Therefore, it is possible to find operators Qi in 
Hi that are similar to the physical operators Oi in Hi, see 
below. 

For the mean-field treatment, one has to find constraints 
which define the Hilbert space Hi in a unique way. As 
shown in Q, this is achieved by means of the operator 
identities 

4o = l-X^W<W=0, (46a) 
r,i 

= $,A,*> (46b) 

- E ®;r,i>kr,i i(I\ia\*>\n = 0- 
r,i,r 
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The representation m, ; r,r' of local operators in Hi is 
readily given by 



mi;r,r> =^<i>\-,r,i<i>i-,r' ,i 



(47) 



This result leads to the representation 

a ;loc = J2 E l°r,r> E $i;r,Ar>,i ( 48 > 



r,r> 



of the local Hamiltonian (f2|i. 

In order to set up the effective Hamiltonian H, one 
needs representations cj of fermionic creation operators. 
As shown in [7;|, a conceivable choice for cj would be 



c f = V a a ' c f , 



(49a) 



where 



" viw-n) ' ; '' :/ - /a ' :, '- / ' 



r,r' 



(49b) 

is a bosonic operator and N is the number of spin-orbital 
states \<r) per site. Evaluated on mean-field level, however, 
expression (|49l does not lead to a reasonable energy func- 
tional since it does not yield the correct results in the uncor- 
rected limit. Note that the situation here is different from 
the single-band model since, there, it was possible to find 
improved expressions for the operator cj a , which are still 

similar to the physical operators cj a . In case of the oper- 
ators (|49l , it seems to be unfeasible to improve them ac- 
cordingly. Instead, Lechermann et al. introduce the follow- 
ing 'improved' expression for d49l . which, though leading 
to reasonable results on the mean-field level, is mathemat- 
ically not similar to cj a , 

€ = EEE ^Un m 

r,r> 7 

*$lriKr'i'Mi-~,,<r> , (51) 



where 

M 



1 



-[A (p) A (h) + A (h) A (p) ]" 1/2 j ( ■ (52) 



and 



r,i,i' 

= E ^rjkrj4l'h*4,M ■ (53b) 



r.i.i' 



Note that 4jw and ^M*,* 

are considered as matrices 
with respect to the indices er, a' whose elements are bosonic 



operators. The inversion [. . .J -1 / 2 and the square root in 
d52i l are defined with respect to this matrix structure. 

The operators d48l ) and (T50b define an effective Hamil- 
tonian 

a E n:f E E C (C ) f 4 Ay + E 



i,j o",er' 7,7' i 

(54) 

which can now be evaluated on mean-field level, i.e., by 
replacing the operators cffip 1 by their corresponding am- 

(*) 

plitudes <p\.pj. These amplitudes then serve as variational 
parameters. 

In Ref. 0, it was shown that the energy functional that 
results from the mean-field treatment of the constraints ( l46l ) 
and of the Hamiltonian ( f54l ) agrees with the Gutzwiller 
variational results introduced in Refs. J4]|5]]. 

5 A new Slave-Boson Theory for Multi-Band Hub- 
bard Models As discussed in the previous section, a gen- 
eralisation of the Kotliar-Ruckenstein scheme for the in- 
vestigation of multi-band models faces significant prob- 
lems, which, up to now, have not been solved satisfactorily. 
Here, we show that, due to the enormous flexibility of the 
slave-boson approach, it is, in fact, relatively easy to repro- 
duce the Gutzwiller energy functional for multi-band Hub- 
bard models. To this end, however, one has to approach the 
problem in a different way than Kotliar and Ruckenstein. 
In the first part of this section, we introduce our new slave- 
boson scheme by reconsidering the single-band model. In 
the second part, we show that our new approach can be 
easily applied to multi-band models including those with 
superconducting ground states. 

5.1 The Single-Band Model For our alternative for- 
mulation of the slave-boson theory, we introduce the oper- 
ators 



where 



and 



1)1: 



it 



n 1 



71=1 



(55a) 



(55b) 



(55c) 



is the projection operator onto the vacuum state of the bo- 
son created by The operator (I55ab therefore projects 
onto the sub-space with exactly one boson in the state \i; I) 
occupied. 

As pointed out before, there is a large amount of arbi- 
trariness in the choice, e.g., of the constraints (f27T >. Instead 
of those equations, we can also work with 

Fifi = 1 - E ™iJ m ¥,i = ' (5 6a ) 
1 



J m f:I 







(56b) 
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where rhi-j has been defined in Eq. ((H). Note that the Hilbert 
space Hi, given by the basis (130) . is already uniquely de- 
fined by the first constraint, equation ( I56ab . In Hi, how- 
ever, the second equation is equally valid, i.e., we have 
F h<J 1/^ = for all states (130) . 

In the Kotliar-Ruckenstein scheme, the operators rhi-i 
are chosen as nf.j, c.f. equation ((28). In our approach, we 
work with 

mi-i = rhi-imfj ■ (57) 

Since we are only interested in the derivation of an ap- 
proximate ground state-energy functional, we avoid func- 
tional integral techniques here and employ the variational 
wave function 

|< B ) = \#o) ® |*b) (58a) 
where |t? ) is a fermionic single-particle product state and 

\Wf) = Y[Di\0) (58b) 

i 

a coherent bosonic state with 
I 

By construction, ( I58bb is a normalised eigenstate of <j> imJ 
with eigenvalues see, e.g., Ref. Ifl9l . In addition, we 
have 

(e»;j>i?B = exp (-|</v| 2 ) , (59a) 
(0i,A,i)*f = |^ ; i| 2 n ex P(-|^'| 2 ) (59b) 
which leads to the expectation value 

where we introduced 

%r = ¥vII exp HW;J'l7 2 )- (61) 
/' 

A comparison with the corresponding result in the Gutz- 
willer theory J4] shows that the Gutzwiller variational pa- 
rameters Xj correspond to the 'renomalised' bosonic am- 
plitudes 

^i;I=Xi;I ■ (62) 

The constraints ((56), evaluated on mean-field level (i.e., 
using the wave-functions d58ab ), have the form 

l = Z)l^/| a <^i;i>»o . ( 63a ) 
i 

n\ a = \$i-„\ 2 {rhi-„)v + |^;i2| 2 (m i;12 }^ (63b) 

which is in agreement with equations (f4TT > if we equate nf.j 
(in the Kotliar-Ruckenstein scheme) with (rhi-i)& 
(in our new slave-boson scheme). 



Finally, we choose the operators cJ. CT in Hi as 

e\,a = %A,c > —i j tr = sJ l0 -Ci, CT ; (64a) 

where 

= *>l.ii<>,,,>',^ + 0\J v3 {l - <ii, 9 ) . (64b) 

Note that here, unlike in the Kotliar-Ruckenstein scheme, 
the operators q^l contain both fermionic and bosonic de- 
grees of freedom. Evaluated with the wave function ( I58al ). 
one finds 

(65) 

which agrees with equations d37) on mean-field level. How- 
ever, the expectation value of a hopping operator is the 
same as in the Gutzwiller theory, 

only if we neglect the 'three-line' contributions 

(^Va^-^;A0) 2 l<^ £ i,ff>*o| 2 (<A«r>*o * . 

(67) 

These terms emerge when the fermionic expectation value 
{hi,gc\ a c- a fij t v)& is evaluated by means of Wick's the- 
orem, 

(Kvc\,<T%Ai,vh = n i,a n< jM4,a d j,a)^o (68) 
M(4,* & j,v)vo\ 2 (4,a d j,a)*o ■ 

Since the three-line terms vanish in the limit of infinite spa- 
tial dimensions, our slave-boson approach yields the same 
variational ground-state energy as the Gutzwiller theory. 

5.2 General Multi-Band Models A generalisation of 
our new slave-boson scheme for multi-band Hubbard mod- 
els is straightforward. We work with an arbitrary set of lo- 
cal multiplet states \r) { , which define the basis ((44) of a 
Hilbert space Hi. For these states, we introduce bosonic 
operators rhf. r and §yl as in ((55) only with I replaced 
by r. 

As a generalisation of ((56), we work with the constraints 

F lfl = 1-1] '»,:/ '"!?/ = , (69a) 
r 

ha,*> = 4,a%a' ~ 4,*%*' ™^™|r = > ( 6 9b) 

r 

which yield an alternative way to define the Hilbert space 
Hi. 

The operators rhi-r r> are properly represented in Hi 

iM'.r.r = mi-,r,r>0i. r ( ) i-,r> ■ ( 7 °) 
An evaluation of these operators on mean-field level, i.e, 
by means of a wave function OH) , with 

Di = Y[cxp ((pi;r4>lr ~ Vl-rk-.r) , ( ? 1) 
r 
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leads to 

mi-r.r' = (mi-,r,r>)&™ = (mi-r,r')^t,r^i;r' (72) 
with 

d t , r = ^r[] cx P(-|^rf/ 2 ) • ( 7 3) 

r> 

A comparison with Eq. <fT9b reveals the correspondence 
of the variational parameters Aj ; r in the Gutzwiller theory 
and the amplitudes fl^r in our new slave-boson mean-field 
approach. 

An evaluation of the constraints (l69l l on mean-field level 
leads to 

1 = ^2K,r^ar(mi;r)»o > ( 74a ) 
r 

r 

which matches the Gutzwiller constraints, equations dT~8T >. 
Finally, we define the operator 

with 

C' = E^^'< r l^l r ')E T ^ T /',r^7,'/' ■ (75b) 
r,r' 

Here, we dropped the lattice-site index i and use the oper- 
ator Hf j, defined in Eq. d23l l. The operator cj a in Hi is 
similar to the physical creation operator cj because the 
sums over er' and /, 7' in equations d75l l are just a compli- 
cated expression for 

\ r ) ( r 'l = EE T i,rTr,r'Hf, r cl . (76) 

cr' 

A mean-field evaluation of d75bl ) leads to the renormal- 
isation matrix 

r,r' 

x E T i,r T i',r' (Hij')vo > 

which matches equation d22l in the Gutzwiller theory. As 
in the single-band model, the expectation values of normal 
and anomalous hopping operators agree with those in the 
Gutzwiller theory 

(c (t) r (t) ) t fb 

only if we neglect the contributions with more than one 
line connecting the sites i and j. This is ensured in the limit 
of infinite spatial dimensions, where both approaches then 
yield the same ground-state energy functional. 



6 Conclusions In summary, we have developed a new 
slave-boson scheme for general multi-band Hubbard mod- 
els which, in the limit of infinite dimensions, reproduces 
the results of the Gutzwiller theory. The main advantage 
of our new approach is its exactness up to the point where 
a mean-field approximation is applied. In addition, it auto- 
matically covers the cases of systems with superconducting 
order parameters. 
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